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^ . Abstract 

Q_ij The dipole moments of the leptons and quarks are matrices in flavor space, 

which can potentially reveal as much about the flavor structure of the theory 
as do the mass matrices. The off-diagonal elements of the dipole matrices lead 
to flavor-changing decays such as fi — ► cy, while the imaginary parts of the 
diagonal elements give rise to electric dipole moments. We analyze the scaling 
of the leptonic dipole moments with the lepton masses in theories beyond the 
standard model. While in many models the dipole moments scale roughly as 
lepton mass, it is shown that simple models exist in which the dipoles scale as 
the cube of the mass or in other ways. An explicit example with cubic scaling is 
presented, which is motivated on independent grounds from large angle neutrino 
oscillation data. Our results have great significance for the observability of the 
electric dipole moments d e , d^, d T , and the rare decays [i — ► ej, and r — ► ^7 
and will be tested in several forthcoming experiments. 



1 Introduction and summary of results 



Most of the parameters of the standard model are masses and mixings of the quarks 
and leptons. These masses and mixings exhibit patterns that suggest to many the- 
orists that there exist underlying "flavor" or "family" symmetries. Over the years 
many models of fermion masses and flavor symmetry have been proposed, some of 
them quite elegant. However, it is impossible without further experimental facts to 
decide among the various theoretical ideas. A determination of neutrino masses and 
mixings will certainly be enormously important for solving the "flavor problem" , but 
probably even more experimental facts will be required. 

Another window on the flavor problem could be provided by the dipole moments 
of the quarks and leptons, since not only the fermion masses but the dipoles as well 
should reflect the underlying flavor symmetries of the theory. In fact, potentially 
there could be as much information about flavor in the dipole moments as in the 
masses and mixings. Given N fermions fi having the same Standard Model quantum 
numbers, their anomalous dipole moments, like their masses, form an N x N complex 
matrix, which we shall denote as D^. This matrix is the coefficient of the effective 
operator f? L Ca tlu F llu f? L . Whereas in the "physical basis" of the fermions the mass 
matrix becomes real and diagonal, the anomalous dipole moment matrix would be 
expected to remain complex and, in general, non-diagonal. 

The imaginary part of the diagonal elements of the anomalous dipole moment 
matrix corresponds to the electric dipole moments of the particles: Im(Dn) = d^. 
The off-diagonal elements of the anomalous dipole moment matrix are the "transition 
dipole moments", which can be either magnetic or electric: Re(Dij) = faj, and 
Im(Dij) = dij. The transition dipole moments lead to decays such as /i — > 

In this paper we shall focus on the dipole moments of the charged leptons. These 
are quite interesting from the experimental point of view as it is likely that there 
will soon be great improvements in the limits on several of these quantities. The 
present limit on the electric dipole moment (edm) of the electron, coming from atomic 
experiments, is 4.3 x 10 -27 ecm ||. A recently approved experiment on PbO expects to 
push this down to the level of 10 -29 ecm in three to five years, and hopes eventually 
to reach even the level of 10~ 31 ecm [[|. The current limit on the muon edm is 
1.1 x 10~ 18 ecm [|J, but there have been recent proposals to push this to 10 -24 ecm 
||. The present limit on the branching ratio of fi — > is 1.2 x 10~ n |J. There 
is a proposal to push this to the level of 10~ 14 at PSI J7J. The forthcoming MECO 
experiment |8j will improve the limit on fx — e conversion in nuclei to the level of 
10~ 16 in branching ratio, which is directly related to the \x — > ej decay rate. The 
Brookhaven experiment E821 is currently in the process of improving the precision 
on muon g — 2 to the level of (few) xlO~ 10 M. The current limit on r — ► fij is 
1.1 x 10~ 6 pi. This limit will also be substantially improved, at least to the 10~ 7 



level, and perhaps to 10 [11]. The current limit on the edm of the tau lepton is 



3.1 x 10 16 ecm fL2| |, which could be improved significantly. 
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In this paper we shall show that there are reasonable and well-motivated models 
in which several of these leptonic dipole moments, specifically D^ e , d^, and D TfJt , 
may be within reach of experiment. If this were so, then the pattern of the leptonic 
edm matrix would be extremely revealing of the underlying flavor structure of the 
theory It is also quite possible, of course, that no lepton edm or transition dipole 
moment is large enough to be seen. Certainly this is the case in the Standard Model, 
where these quantities are highly suppressed. However, if there is non-trivial flavor 
structure at scales between a few hundred GeV and a few hundred TeV, then some 
of these moments may be observable. If there is no new physics all the way to 
about a thousand TeV, then there is very little chance for seeing these edm's or the 
rare leptonic decays in the forthcoming experiments. However, new physics such as 
supersymmetry is anticipated not much above a few TeV, so chances for seeing these 
effects are great, as we will argue. 

In many models there is a close relationship between the anomalous dipole moment 
matrices and the mass matrices of the fermions. Both kinds of matrix connect the left- 
handed fermions to the right-handed fermions. Moreover, at the diagrammatic level 
the Feynman diagrams that contribute to dipole moments and masses are related. In 
particular, removing the external photon from a dipole moment diagram produces a 
mass diagram, and adding a photon to a mass diagram of one or more loops produces 
a dipole diagram. In fact, in many models the anomalous dipole moments of the 
quarks and leptons are roughly proportional to their masses. 

In the minimal supersymmetric Standard Model (MSSM), for instance, there are 
well-known one-loop contributions to lepton and quark edms coming from gaugino 
loops [nj. In those diagrams a lepton (or quark) emits a virtual gaugino to become 
a virtual slepton (or squark) and then reabsorbs the gaugino, coupling in the process 
to an external photon line. Since a dipole moment operator involves a chirality flip of 
the fermion, there must be an insertion either of the mass of the external fermion or 
of the "left-right" mass of the corresponding virtual sfermion. However, the left-right 
sfermion masses, m 2 LR , are proportional to the masses of the corresponding fermions 
in the MSSM. Thus the edms of the quarks and leptons come out being approximately 
proportional to their masses. 

Proportionality of edms and masses also occurs in multi-Higgs models. In models 
with more than one Higgs doublet there can be CP violation in the Higgs propagators 



that can lead to measurable edms ||14|| . The dominant contribution to the edms of 
the lighter quarks and leptons come from two-loop diagrams that involve one power 
of the Higgs Yukawa couplings [[ilj . These contributions are therefore approximately 
proportional to the mass of the quark or lepton. (There are also one-loop diagrams, 
but these are cubic in the Higgs Yukawa couplings and are therefore negligible for 
the lighter families. These one-loop contributions may however become significant 
for the muon edm for some range of parameters, especially large tan/5 ||16|| .) 

This feature that the anomalous dipole moments scale roughly as the first power 
of the fermion mass is characteristic of many models. If it does in fact hold in Nature, 
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then it has several important implications. First, there would be a limit on the muon 
edm coming from the present limit on the electron edm: ~ (m^/m^de < 8.8 x 10~ 25 
ecm. This is at the margin of what could be seen by the proposed experiments. 
Second, one would conclude that d T and D Tfl are unobservably small. For example, 
one would have d T ~ (m T /m e )d e < 1.5 x 10~ 23 ecm, which is 0.5 x 1CT 7 of the present 
limit and is well beyond the reach of currently planned experiments. 

However, the anomalous dipole moments of the leptons and quarks do not neces- 
sarily scale as the first power of the masses. In this paper we shall describe models in 
which they scale as other powers, in particular as the cube or the square. The idea 
that electric dipole moments might scale as the cube of the mass is not new. As al- 
ready noted, the one-loop contributions to the fermion edms in multi-Higgs models do 
tend to scale with the cube of the mass, and these were for a time naturally thought 
to be the leading contributions. However, it was then pointed out Jl5[ that two-loop 
diagrams which scale as the first power of mass dominate for the lighter families. 

We present an explicit model in Sec. 2 and 3 where the dipoles obey a cubic scaling. 
The model is motivated on independent grounds, viz., an explanation of solar and 
atmospheric neutrino oscillation data in terms of large angle neutrino oscillations. 
In this model we find that the anomalous dipole moment matrix scales roughly as 
follows: 



d T 



( (m e /m T f + cf-{m e /m T ) 



(m e /m T ) 2 



\ 



\ {m e /m T y (m M m T ) 1 J 

Here the dashes indicate elements that are negligibly small. One sees that the diagonal 
elements — and therefore the edms of the particles — scale as the mass cubed, while 
those elements below the diagonal are suppressed relative to the diagonal elements 
by factors that go as the mass ratios (coming from mixing angles). 

In Eq. (1), the dipole moment of the electron has a term that is denoted as 
c^{m e /m T )d T . Since the electron mass is so small, one has to take into account two 
loop diagrams for the dipole moments which may have a linear scaling with the lepton 
mass Hl5| . Typically, we find that the coefficient c is of order unity. With c = 0(1), the 
two loop diagram with linear scaling will win over the one-loop diagram with cubic 
scaling. But this is significant only for the electron dipole moment. It is also possible 
that these two-loop diagrams are suppressed in some models, in which case c ~ 0, 
although in the explicit models that we have constructed, this turns out not to be 
the case. 

Consider the case of c ~ (i.e., the two-loop diagrams being suppressed). The 
resulting pattern is interesting from the point of view of several kinds of experiments. 
If we suppose that the transition dipole D^ e is just at the current limit coming from 
B(fi — > e'y) < 1.2 x 10 -11 , namely \D^ e \ = 4.5 x 10~ 27 ecm, then the pattern in Eq. 
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(1) with c = would imply the following values of other moments: d e ~ 10 31 ecm, 

5 x 10~ 21 ecm, and \D Tfl \ ~ 3 x 10~ 22 ecm (which corresponds to 



dp ~ 



10 



~ 24 ecm, d T 
/I'j) ~ 10' 



B(t — > /I'j) ~ 10 ). It should be emphasized that these are very rough estimates, 
which could well be off by an order of magnitude in the dipole moment or two orders 
of magnitude in a rare decay rate. Making allowances for this, one can say that if 
the pattern in Eq. (1) is correct then: (a) d e is probably too small to be seen in 
the near future, though experiments on molecules like PbO have a chance of seeing it 
eventually. (And even if not seen, the failure to see it together with the observation of 
would confirm a more-than-linear scaling of the edms with mass.) (b) d^ should be 
at or near the level which recently proposed experiments can reach. And (c) r — > 717 
should be very close to the present limit. 

In the more realistic case of c = 0(1), the above estimates will still hold, except 
for the electron edm. d e should be also observable, at the level of 1 x 10~ 27 ecm, if 
the decay fi — > ej is near the present experimental limit. 

In section 4, we shall describe, in less detail, a model that illustrates quadratic 
scaling of the anomalous dipole moments with mass. In that model one finds 



d T 



I {m e /m T ) 2 (m e mp/m 2 ) {m e /m T ) \ 
(mjrip/ml) (m^/m T ) 2 (m M /m T ) 



(2) 



\ KK) {m^/m T ) 1 / 

In this case the two-loop contribution to d e which would scale linearly with the 
electron mass is at best of the same order numerically as the one-loop correction and 
is not shown in Eq. (2). If /i — > is saturated at the present limit, this pattern 
would lead to the following predictions: (a) d e ~ 2 x 10~ 29 ecm, (b) rf^lx 10~ 24 
e-cm, (c) d T ~3x 10~ 22 ecm, (d) \D Tfl \ ~2x 10~ 23 ecm (corresponding to B(t — > 
/j^y) ~ 4 x 10~ 7 ). Again, we see that these predictions can be well tested directly in 
the near future. 



2 Flavor symmetry and the dipole moments 

Before describing a specific model that leads to the pattern shown in Eq. (1), we will 
describe in more general terms how different kinds of scaling of the lepton edms can 
arise in simple models based on abelian flavor symmetry. 

Many models of fermion masses explain the interfamily mass hierarchy by appeal- 
ing to spontaneously broken flavor symmetry. The idea is that certain elements of 
the mass matrix would vanish if some abelian symmetry or symmetries were exact, 
and are therefore suppressed by powers of the expectation values that violate those 
symmetries. The fields whose expectation values violate flavor symmetry are often 
called "flavons" . If F denotes a flavon field, then different elements of the mass ma- 
trix could be suppressed by different powers of e = (F)/Mp, where Mp is a scale 
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characterizing the underlying flavor physics. The absolute scale of Mp and (F) could 
be anything from just above the weak scale up to the Planck scale. All that matters 
for explaining the fermion mass hierarchy is that their ratio e be smaller than, but 
of order, unity. Some models have several flavon fields, and have the fermion masses 
depending on several small parameters £j. 

In a typical model of the kind we have been describing the fermion masses arise 
from tree-level diagrams similar to that shown in Fig. 1(a) [fL7|| . One of the boson 
lines attached to the fermion line is a Higgs doublet whose vacuum expectation value 
breaks SU(2)l x U(1)y- The others are flavon fields, whose vacuum expectation 
values give rise to the requisite powers of the small parameter (s) £j. It is evident that 
if two of the boson lines are tied together to make a loop, and an external photon line 
is attached, the resulting diagram (Fig. 1(b)) gives a non-zero dipole moment to the 
fermion. If the flavor scale is near the unification scale or Planck scale, the resulting 
dipole moment is negligible. However, if the flavor scale is near the weak scale, a 
significant dipole moment can result. One sees, moreover, that the dipole matrix will 
have a close relationship to the mass matrix. 

Consider the following simple toy model as an illustration. This model will be gen- 
eralized in the next section to a fully realistic model that also explains large neutrino 
mixing angles as suggested by recent data from solar and atmospheric neutrino exper- 
iments. The toy model is over-simplified in that it has vanishing inter-generational 
mixings. Suppose one has the ordinary three families of leptons (£i) and three families 
also of singly- charged vectorlike leptons (Xi): 

= Li, £ iL , X iL , X iL . (3) 
Suppose the following Yukawa couplings and mass terms: 

C Yuk = £ h{Utt L )rt + Y: ff{X- L X^ L )h K + £ ffiXrJ^S,. (4) 

i i,K i 

This structure can be enforced by a simple flavor symmetry. The field H is the usual 
Higgs doublet, while and hx are singlets. Call fi(H°) = rhi, Ex//^ {Hk) = (F), 
and fi(Si) = Mi = e^ 1 (F i ), where it is assumed that < 1 and <C Mj, (F,). 
Then there is a mass matrix of the form 

^ X * l) {m\ (F i )){x+)- (5) 

This matrix has three large eigenvalues that are of order Mj, and three small eigen- 
values that are of order rhi. The three light eigenstates are to be identified as the 
physical e, /i and r. Neglecting terms of order (rhi /Mi), the large entries in the 
mass matrix are diagonalized by a change of basis of the positively charged leptons: 
I'll = sin9i£f L — cos6iXf L , and Xf[ = cos Oi£f L + sin^X^, where tan#« = e^. (We 
have implicitly chosen a phase convention wherein the entries in the mass matrices 
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have all been made real. The couplings of the leptons to the scalar fields will however 
remain complex in this basis.) This would give effective mass terms for the light three 
fermions of the form mi(£~ L £fl), with 

(Ft) . . 

m = -j^m = rriiei, (6) 

where we have used sin#j = tan#j = 6j. We suppose that the parameters rhi are all 
of the same order and that the fermion mass hierarchy comes from a hierarchy in the 
ef e 1 < e 2 < e 3 . 

Diagrammatically, one has at lowest order in the the tree-level diagram in Fig. 
2(a). This would give an effective mass term for the light fermions of (fi{H))(ff {Si})^ 1 
CE K ff{h K )) = m i M~ 1 {F i ) = rhiei, which of course agrees with the result of diago- 
nalizing the matrix. 

If an external photon line is added to the diagram in Fig. 2(a) and the H and 
hx boson lines are tied together through the quartic term \ K H^ Hh* K h K , which must 
certainly be allowed to exist by symmetry the diagram shown in Fig. 2(b) results. 
The contribution to the dipole moment of the i th light lepton [£ x = e, £ 2 = A*, and 
£3 = r) will be of the form 

d i = j^E( eX Kfdn(H)(h K )MJ(m 2 H ,m 2 hK ,M?), (7) 

where I is the function with dimensions of M~ 4 resulting from doing the momentum 
integral, given by I(r 1 M,r 2 M,M) = |M~ 4 (ri - r 2 ) _1 ([l - r\(l +lnr?)]/(l - n) 3 - 
(ri — > r 2))- Let us suppose that the largest mass in the loop is that of the heavy 
virtual fermion (which Eq. (5) shows to be approximately Mj). Then the integral 
is given by / = (M 4 )- 4 (| + \{m\ + m^)/Mf)), Recalling that Z K f t K (h K ) = 
define J^xXxf^ihR) = Aj(Fj). If one neglects the K dependence in the momentum 
integral J, the dipole moment of the i th light lepton can be written as 

It is technically natural to assume that the CP-violating phase of A, is of order one. 
It should be noted that there will be such a phase in the physical basis of the leptons, 
in general, if there is more than one type of boson hx- The point is that the mass 
of the i th light lepton, which is a real number in the physical basis, is proportional to 
^Kfl^ihK), whereas the dipole of that lepton depends on SA^-//^ (hx), which has no 
reason to be real in the same basis. Thus the electric dipole moment of the i th lepton 
is of the same order as the result in Eq. (8) 

One can see that this edm can be quite large for reasonable values of the pa- 
rameters. For example, if M 2 = 3 TeV and X 2 — 1, then the edm of the muon 
would be of order 3 x 10~ 24 ecm. The scaling of the edms with mass depends 
upon how the family mass hierarchy is assumed to arise. This hierarchy is given 
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by mi/rrij ~ ei/ej = {(Fi) / (Fj))(Mi/Mj)~ 1 . Thus, it can come from (a) a hierarchy 
in the {Fi), with all the M, being comparable, (b) a hierarchy in the Mi, with all the 
{Fi) being comparable, or (c) a combination of these. If (a) is the case, then Eq. (8) 
would imply that the edms of the leptons scale approximately linearly in lepton mass. 
However, if (b) is the case, then the edms scale with the cube of the masses, since 
mi/rrij ~ Mj/Mi and thus di/dj ~ (mi/mj)(Mj/Mi) 2 ~ (m^/m.,) 3 . 

Since the one-loop dipoles scale as the cube of the lepton masses if the hierarchy 
is in Mi, one should examine whther two-loop diagrams induce larger contributions, 
especially for d e . We find that this is indeed the case, as in multi-Higgs models. 
Diagonalizing the mass matrix of Eq. (5) to lowest order in e, and m^/M, will lead 
to the following Yukawa couplings of the light fermions to the scalars: 



£ Yuk = iT^ifiH + i- L ifl^f t K h k + H.C. (9) 

In the physical basis, e^mj are real, but f± are complex in general. Thus, an edm of 
the electron can be induced through a two-loop diagram with Hk propagating. The 
dominant contribution to d e will arise when hx decays into a real and a virtual photon. 
This happens, for example, by a diagram in which Hk emits an X^ L X^ L pair which 
closes by emitting two photons, or one where hx becomes an H (through quartic 
scalar couplings), which then decays to two photons through a top quark (and other) 
loop. All these diagrams, while they are suppressed by an additional loop factor, 
scale linearly with the lepton mass and can be dominant over the one-loop diagram. 
For rrih K — 300 GeV, d e turns out to be of order 4 x 10~ 27 ecm |Tj|, if the relevant 
phases are all of order one. 



3 A realistic model with mixing of charged lepton 
families 

The simple example just discussed is perhaps too simple in that it does not have 
any mixing among the families. It is true that the leptonic mixing seen in neutrino 
oscillations might arise from the neutrino sector rather than the charged lepton sector. 
However, the CKM mixing of the quarks and the likelihood that quarks and charged 
leptons are unified in some way suggest otherwise. In this section we present a 
unifiable model of lepton mass in which the charged lepton mass matrix is off-diagonal 
in the flavor basis. This model is actually a realization of an idea proposed several 
years ago to explain both the fermion mass hierarchy and the largeness of neutrino 
mixing angles ||18|| . The model is thus well- motivated on grounds having nothing to 
do with dipole moments. But it naturally leads to the pattern given in Eq. (1). 
The model has the following leptons and scalar fields: 
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field: t~ if X~ X+ H Si h K 

F charge : $ -ft -q t ^ 

where F is an abelian flavor symmetry. This symmetry allows the following Yukawa 
terms for the leptons: 

c Yuk = £ Mere})H° + £ f*{XrX+)h K + £ fS(Xre})Si. (11) 

The crucial difference between this model and the toy model discussed in the last 
section is that here the Yukawa coupling is a matrix. The coupling f$ is also a 
matrix, but this is of less significance, as we shall see. It will simplify the discussion 
without affecting the qualitative conclusions if we assume that f$ is diagonal: f$ = 
ffSij. This assumption will be relaxed later. 

Let us define the matrix ra^- = fij(H°). There will be a 6 x 6 mass matrix of the 
charged leptons, which has the same form as in the toy model (cf. Eq. (5)) except 
that rhij is here non-diagonal: 

i^x- L ){ tl. )( v't )■ (12) 



(F;}6 



■ij 




The quantities Mj and (Fj) appearing here are defined in the same way as in the toy 
model, and = (Fj)/Mj. After diagonalizing the large, i.e. 0{M i ) ) elements in the 
mass matrix, the light leptons (namely £~ and = X^ — e^) have an effective 
mass matrix given by 

(F) 

m ij = ™*3~m7 =™v e r ( 13 ) 

As in the toy model of the last section, we assume that all the f\j and hence all the 
ffiij are of the same order, and that the fermion mass hierarchy comes from the €f 
€i e 2 e 3 . From Eq. (13) it is obvious that the 3x3 light lepton mass matrix has a 
column form. That is, the first, second and third columns of the matrix have elements 
that are, respectively, small, medium, and large. That means that in diagonalizing 
it, the rotations that are done from the left (i.e. of the £~ L ) are of order one, whereas 
those done from the right (i.e. of the ifl) are suppressed by the hierarchy factors 
e j/ e k ~ fnj/mki j < k. This means that there would be large leptonic mixing angles, 
as is observed in atmospheric neutrino oscillations, and perhaps also in solar neutrino 
oscillations. If one were to embed this model in SU (5) then the mass matrix of the 
down quarks d, s, and b would have a similar structure, except transposed (since 
SU(5) relates M L to M%, where M L and M D are the mass matrices of the charged 
leptons and down quarks). As the down quark mass matrix would thus have a row 
structure instead of a column structure, it would be the mixing angles of the dn that 
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would be of order rrii/mj, while those of the d\ L would be of order one. Thus, such a 
model would explain the fact that at least certain of the leptonic (MNS) mixings are 



large while the quark (CKM) mixings are small 

Returning to the leptons, it is easy to see that their dipole matrix also has a 
column form. In fact, the same calculation that led to Eq. (8) gives 

Both the dipole matrix and the mass matrix of the light leptons have a column form, 
and the corresponding columns of the two matrices are in fact proportional to each 
other. However, the hierarchies among the columns are different in the two matrices. 
(Compare Eqs. (13) and (14).) The case we are most interested in is where the 
hierarchy in the e« parameters comes from a hierarchy in the Mj, with the (Fi) all 
being of the same order. Then the columns of m^- are in the ratio e\ : e<i : 
while the columns of are roughly in the ratios ef : e| : e|. This shows that 
d e : dfj, : d T = rn? e : : rr? T . 

To find the magnitudes of the off-diagonal elements of is simple. The first 
thing to be noticed is that because of the fact that the columns of m^- and are 
proportional to each other, the same unitary transformation acting from the left will 
make zero all the elements above the diagonal in both matrices. This is a crucial fact, 
since otherwise there would be large off-diagonal elements in Dij. For instance, in the 
matrix in Eq. (14) the (12) element is of the same order as the (22) element, which 
leads to the danger that d^ ~ |Z) eM |. That would mean that the limit on fi — ► c-f 
would put a severe constraint on the muon edm. As it is, the elements above the 
diagonal can be made zero simultaneously in both matrices, 

To complete the diagonalization of m^- requires a unitary transformation from the 
right. Because there is a strong hierarchy among the columns of rriij this unitary 
transformation involves rotations C/y, i < j, that are of order e$/ej <C 1. This will not 
eliminate the elements of D^ that are below the diagonal, since the hierarchy among 
the columns of is different. Rather, the result of this unitary transformation will 
be to make Dj iy i < j, be of order (ei/e^Djj. This gives the form shown in Eq. (1). 

An interesting prediction that arises from the nearly triangular form in Eq. (1), 
aside from the ones already discussed in the Introduction, is that in the decay fi~ — > 
e~7 the electron will be almost purely right-handed. A similar remark applies for the 
decay t~ — > 

Up to this point we have oversimplified the analysis of this model by assuming 
that the Yukawa matrix f$ appearing in Eq. (11) is diagonal. The symmetries do not 
require it to be diagonal, so this is an unnatural assumption. However, eliminating 
this assumption makes no qualitative difference to the form of the dipole matrix. 
Suppose that f$ is allowed to take a general form. Then the 6x6 mass matrix of 
the leptons is 
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( 



fhn 



Mij (F z )5,, 



) 



( 




) 



(15) 



where My = ffj(Si). As before, we will assume that the hierarchy comes from 
My = ffj(Si) rather than from (Fj). Since f$ is assumed to be an arbitrary matrix, 
this means that the hierarchy comes from (Si). In other words, (Si) 3> (S 2 ) 3> {S3). 
This means that has a row structure, with the first row being much larger than 
the second, which in turn is much larger than the third. By rotating among the £^ L 
one can go to a basis where is triangular, with zeros above the diagonal. This will 
also change the form of the Yukawa matrix f^, but since that matrix is not assumed 
to have any special form anyway, this makes no difference. Thus without loss of 
generality, one can take to have this triangular form. To completely diagonalize 
it then requires a rotation among the X~ by angles that go as (Sj) / (Si) ~ Cj/ej, % < j ■ 
This rotation will have the effect of making the coupling of X~ to Xj~ non-diagonal. 
This is the only change from our previous analysis. However, the off-diagonal elements 
thus introduced are suppressed by the small hierarchy factors e^/e^, i < j. One can 
easily trace through the effect of these small off-diagonal elements, and one finds that 
they give rise to a mixing of the i th and j th columns of the dipole matrix that is of 
order e«/ej, i < j. This does not change the form of the dipole matrix given in Eq. 



As in the previous section, there are two-loop diagrams that are suppressed by 
additional factors of (a/47r) ~ 10~ 3 , but these are still more significant to d e . The 
dominant diagrams involve the exchange of hx fields. The coupling matrix of hx to 
the light fermions are neither diagonal nor real in the physical basis for the leptons. 
As in the toy model, we would expect d e ~ 4 x 10~ 27 ecm if the masses of hx are 
of order 300 GeV. These two-loop diagrams have negligible effects for the muon and 
the tau lepton. 

We have a model of mixing that is satisfactory for the leptons. As already men- 
tioned, it is possible to embed this structure in an SU(5) model. The simplest way 
to do this would be to add to the usual three families of 10; + 5; three vectorlike 
pairs 10^ + 10^. The Higgs Si and h K would be SU(5) singlets, while H would be 
a 5. The charge assignments in Eq. (10) and the couplings in Eq. (11) would be 
generalized in the obvious way. This would lead to a mass matrix for the down quarks 
that was simply the transpose of that for the charged leptons, just as in the mini- 
mal SU (5) model. A realistic SU (5) version of this model would have to introduce 
non-minimal features, such as larger Higgs multiplets or higher-dimensional effective 
operators contributing to light fermion masses. That would allow Clebsch coefficients 
to appear in the mass matrices which differed for the down quarks and the charged 
leptons. Nonetheless, one would expect that the down quark mass matrix, Mp, would 
be closely related to the transpose of the charged lepton mass matrix, Mj, even if 
not exactly equal to it. As noted above, this would elegantly explain the curious fact 
that all the CKM quark mixing angles have been observed to be small, while at least 



(1)- 
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some of the MNS leptonic angles are large. For M L has a column structure, while 
Md has a row structure. 

A question arises, however, whether the model can be extended to include quarks 
without running into a problem with excessive flavor changing in the quark sector. 
In particular, will one-loop diagrams involving the scalar fields H, Si, and hx lead to 
excessive mixing? The answer turns out to be no. Consider, for example, 

an SU (5) version of this model where one takes the parameters to have the following 
orders of magnitude: / - ~ KT 2 , fg ~ 1, f§ ~ 1, (H) ~ 300 GeV, (h K ) ~ m hK ~ 300 
GeV — 1 TeV, (S 3 ) ~ 300 GeV — 1 TeV, (S 2 ) ~ (e 3 /e 2 )(S 3 >, and (ft) ~ (e 3 /e 1 )(S 3 ), 
where e\ : e 2 : e 3 ~ m e : m M : m r ~ : m s : m^. It can be shown that none of the 
box diagrams involving the scalar fields H, h K , and Si, and the new heavy vectorlike 
fermions contribute excessively to either the real or imaginary parts of M 12 in the 
neutral kaon system. Some of the diagrams are suppressed by small mixings of order 
md/rrib and rris/mb or powers thereof, or by the large masses of the vectorlike fermions, 
especially Xf and X 2 ~ ', or by the small values of the fy, or by combinations of all these 
factors. This fact, namely that one can have in this model one-loop dipole moments 
that are of interesting magnitude — including flavor- changing ones — without having 
excessive one-loop flavor changing in the kaon system, is non-trivial. 

It might be suspected that the mixing of the light fermions with the heavy vector 
fermions could lead to large flavor-changing Z couplings that can arise at the tree 
level. However, this turns out to be not the case in this model. For example, in 
the lepton sector, there is no flavor changing Z coupling at the tree level in the 
right-handed sector. There are such couplings in the left-handed sector, but they are 
suppressed by (rh/M) 2 in the amplitude, which is at most of order 10~ 6 for ~ 10~ 2 
and M of order TeV. The resulting branching ratio for Z — > t/i decay for example 
is at the level of 10~ 12 , which is too small to be observable, but is not much below 
the current limit for the decay fi — > 3e. 

It should be observed that the 577(5) version of this model described above could 
not be supersymmetrized if we wish to maintain perturbative unification at a scale 
of 10 16 GeV. Having three vectorlike 10 + 10 pairs along with supersymmetry would 
destroy the perturbativity of the gauge couplings around 10 9 GeV. 

4 Other scalings of the lepton dipole matrix 

We have seen examples of linear and cubic scaling of the electric dipole moments 
with fermion mass. It is not difficult to construct models in which the scaling is 
quadratic. For example, consider a model in which there are the following leptons 
and Higgs fields. Leptons: L t = {v iL ,£~ L ) T , £f L , n aL , n c aL . Higgs: H = (H°,H-) T , 
h a = {h a , h~) T , s~. Here n a are neutral leptons which are taken to be Dirac particles. 
It is assumed that the fields h a and s a have vanishing VEVs. Let there be the following 
couplings: H a m a (n c aL n aL ) + H a)i f ia h\(Lin a ) + H ayi f' ia s~(£ + n a ) + E a X a h a H(s~)* . This 
form can be enforced easily with abelian flavor symmetries that distinguish fields with 
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different values of the index 'a'. 

Suppose that the masses of the known leptons £ arise predominantly from the 
loops shown in Fig. 3. Let the heaviest particles in the loops be the bosons h a and s~, 
which are assumed to have masses ~ M a that are in a hierarchy: M3 <C M2 <C M\. 
The largest contribution to the lepton mass matrix m^- therefore comes from the 
diagram in Fig. 3 having a = 3. This gives a contribution to m^- that is of rank 1. 
This would generate the mass of the r lepton. The second largest contribution comes 
from the diagram with a = 2 and also is rank 1. This would give mass to the muon. 



The diagram with a = 1, finally, gives the electron mass [20]. If we assume that the 



only hierarchy in the underlying parameters is in the M a , then it is apparent from 
simple power counting that the lepton masses have a hierarchy that goes roughly as 
m e : m M : m T ~ Aff 2 : M^ 2 '■ M% 2 , since the loop diagram of Fig. 3 is quadratically 
convergent. 

If one adds an external photon line to the loop, it must couple to the boson line, 
since the virtual lepton n a is neutral. The resulting diagram for the anomalous dipole 
moment is quartically convergent. Thus one obtains d e : : d T ~ Aff : M^ 4 '■ 
Mg -4 ~ m 2 : m 2 : m 2 . In fact, it is not difficult to show that the whole anomalous 
dipole matrix has the scaling given in Eq. (2). 

One difference from the model of Sec. 2 and 3 is that numerically the magnitudes of 
the dipole moments are enhanced here by a loop factor, which is roughly (1/167T 2 ) -1 ~ 
100. The observability of the edm is thus greater in this case. 

It is also possible to generate the dipole matrix of Eq. (2) by introducing singlet 
charged leptons and neutral scalars instead of neutral leptons and charged scalars. 

There is an obvious generalization of this model to the quark sector. We can 
introduce singlet quarks with electric charges +2/3 and —1/3. The ordinary quarks 
can couple to these vector quarks through the higgs doublets h a = (h a ,h~) and 
neutral higgs singlets s a which do not acquire VEVs. One-loop diagrams analogous 
to Fig. 3 will induce quark masses as well as dipole moments. If the masses of the 
h a fields (and/or the s a fields) are hierarchical, the observed fermion mass hierarchy 
can be reproduced, with all Yukawa couplings being of the same order. This might 
be an interesting radiative mass hierarchy model for all the fermions, except for the 
top quark. 

In the supersymmetric standard model, if there is non-trivial flavor structure 
in the trilinear A terms, analogous mass hierarchy can be induced, with interesting 
dipole moments [^TJ. In this case, the diagrams will involve exchange of the gaugino 
and squarks. 



5 Conclusions 

The magnetic and electric dipole moment matrices of quarks and leptons can poten- 
tially carry a wealth of information about the underlying flavor structure. In a large 
class of models we have shown that these are experimentally accessible, through the 
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electric dipole moments of the electron, muon and the tau lepton, and through rare 
decays such as \i — > and r — > n^. We have addressed the scaling of the leptonic 
dipole moments with the lepton mass. While linear scaling with the mass is quite 
generic, in a class of models motivated on independent grounds by large angle neu- 
trino oscillation solution to the solar and atmospheric neutrino data, we have found 
interesting scalings, where the dipoles go as the cube or the square of the lepton mass. 

Dipole moments of quarks and leptons have been analyzed in the literature in 
various other contexts. Their scaling with mass in grand unified supersymmetric 



models is typically linear [^, |23| , in contrast to the cubic scaling that we have found. 
It is also possible that significant contributions to leptonic dipole moments arise 
through SUSY exchange from interactions that are responsible for heavy right-handed 



neutrino masses [33]. Such interactions do not lead to any simple power law scaling. 
Experiments in the near future will be capable of telling these theories apart. Indeed, 
very important information about the flavor structure will be gained by the proposed 
improvements in the dipole moment measurements as well as in the radiative decays 
of (jl and r. 
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(H) (F t ) (Fj) 

Fig. 1(a) 

A typical tree diagram for light fermion mass through flavon fields. 




Fig. 1(b) 

Dipole moment of the fermion induced by folding the legs of Fig. 1(a). 
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Diagram that induces lepton masses in the model of Sec. 2. 
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Fig. 2(b) 

Lepton dipole moment contribution in the model of Sec. 2. 
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Fig. 3 

Diagram that induces radiative masses for the leptons in the model of Sec. 4. 
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